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Using a new analytic quantum mechanical method based on Slater's Xa method, we show that a 
fairly accurate estimate of the total energy of a molecule can be obtained from the exact energies of 
its constituent atoms. The mean absolute error in the total energies thus determined for the G2 set 
of 56 molecules is about 16 kcal/mol, comparable to or better than some popular pure and hybrid 
density functional models. 
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PACS numbers: 

The total electronic energy is the central quantity in 
quantum mechanical methods that compute the ground 
state properties of many-electron systems. These meth- 
ods rely on the variational principle and can be broadly 
classified in two categories. The first category is formed 
by ab initio quantum chemical methods which permit sys- 
tematic improvement in accuracy and analytic computa- 
tion of quantum mechanical matrix elements to machine 
precision^ The high computational cost of these meth- 
ods restricts their application to those systems with a few 
atoms. Density-functional-based computational models 
form the second category. These models are based on 
the Kohn-Sham 2 formulation of density functional the- 
ory (DFT) which offers an alternative way to determine 
ground state properties. Today DFT based models have 
become the most popular choice for calculating these 
properties as they provide results of sufficient accuracy 
at reduced computational costs. Unlike ab initio quan- 
tum chemical methods, however, the computational im- 
plementation of almost all DFT models require the use of 
numerical grids- making calculation of matrix elements 
at machine precision practically impossible. Round-off 
error, which grows as the square root of the number of 
points, is eliminated by analytic methods. 4 

The desirable attribute of analytic computation, which 
is accurate to machine precision, is also possible within 
DFT.'' This approach is based on variational and robust 
fits to the orbitals and the effective one-body Kohn-Sham 
potential using Gaussian basis setsA This procedure does 
not alter the density- functional one-particle equations, 
but makes its analytic solution possible. The model has 
its roots in the early density-functional model of Slater, 7 
wherein the exchange potential v x is proportional to the 
one-third power of the electron density p: 
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where a is called the Slater exchange parameter. We 
extended this analytic density-functional model so that 
an atom-dependent exchange parameter a, which scales 
the exchange potential differently for each element, can 



be used in heteroatomic molecules or solids^ Early nu- 
merical attempts to use atom dependent a values had a 
mathematically undefined total energy that results from 
the discontinuity in the potential at the boundary of 
the muffin-tin sphere that enclosed atoms or ions, in the 
Multiple-Scattering Xa method. 9 In our method the to- 
tal energy is both well defined and stationary with re- 
spect to variations of all linear combination of atomic- 
orbital (LCAO) and Kohn-Sham-potential fitting coeffi- 
cients. The energy is a function of basis sets and Slater's 
a for each element, 

E({ai}) = ^<^|/i|0 i >+2(p||p}-(p||p) 
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Here, /i is the one-electron part of the Hamiltonian, p a 
is the spin density and an overbar represents an LCAO 
fit, and gs, the partitioned 3/4 power of the exchange 
energy density: 
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where Dfj [f) is the diagonal part of the spin density ma- 
trix, and the function, 
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contains on, the a within the muffin-tin in the Xa 
method, for the atom on which the atomic orbital i is 
centered. This unique expression for the total energy 
is obtained by adding different fitted expressions for the 
same components of the energy so as to cancel all first- 
order errors in the energy due to all fits. All orbitals and 
fits are determined through variation of Eq. ^ 
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TABLE I: The optimal a values that yield the exact atomic 
energies in the highest symmetry for which the solutions have 
integral occupation numbers for the analytic DFT calcula- 
tions. The exact atomic energies given in the last column 
are from Refill The basis sets are I: 6311G**/RI-J, II: 
DZVP/A2. 

Basis I Basis II Numerical E (a.u.) 

TE 0.77739 0.78124 0.77679 -0.500 

Li 0.79169 0.79211 0.79118 -7.478 

Be 0.79574 0.79614 0.79526 -14.667 
B 0.78675 0.78677 0.78744 -24.654 

C 0.77677 0.77665 0.77657 -37.845 

N 0.76747 0.76726 0.76654 -54.590 

O 0.76500 0.76448 0.76454 -75.067 

F 0.76066 0.76001 0.75954 -99.731 

Na 0.75204 0.75287 0.75110 -162.260 
Mg 0.74994 0.75120 0.74942 -200.060 
Al 0.74822 0.74869 0.74797 -242.370 
Si 0.74539 0.74602 0.74521 -289.370 

P 0.74324 0.74397 0.74309 -341.270 

S 0.74262 0.74350 0.74270 -398.140 

CI 0.74197 0.74272 0.74183 -460.200 



This method inherits from Slater's Xa method the 
physically appealing advantage, which it had over all 
other quantum-chemical methods, that atoms dissoci- 
ate correctly. For these molecules and this range of a 
values all dissociated atoms turn out to be neutral in 
spin-polarized (high-spin), broken-symmetry (the high- 
est symmetry that gives integral occupation numbers in 
fractional-occupation-number) calculations. Using Eq. 1 
and atomic a values is a single approximation, like the lo- 
cal density approximation is a single approximations and 
not a different parameterization for each range of densi- 
ties. In particular, we determine a set of atomic a values 
that yield the exact atomic (EA) energies and use them 
to compute total energies of the G2 sekii of 56 molecules. 

Our calculations use Gaussian basis sets to fit both the 
orbitals and the Kohn-Sham potential. Here, we choose 
two different combinations of basis sets. For orbitals, we 
use the valence triple-C (TZ) 6-311G** basisi£ and the 
DGauss valence double-C basisii set (DZVP2). The s- 
type fitting bases are obtained by scaling the exponents 
of the s part of the orbital basis, by two to fit p, by 2/3 
to fit p 3 j and by 4/3 to fit p 3 ■ For the non-zero angular- 
momentum components the resolution-of-the-identity-J 
(RI-J)A£ and A2i4 basis sets are used without scaling 
to fit the Coulomb potential as well as the exchange- 
correlation part of the Kohn-Sham potential. Thus, two 
combination sets 6-311G**/RI-J and DZVP/A2 of bases 
were used for obtaining the atomic a values. The a op- 
timization was performed using the PERL scripts. 

We determine of a values that give the exact total 
electronic energy for each atom in the molecule. That 
set of a values is given in Table [I] for each basis set. 
The "exact" atomic energies^ that these values of a re- 
produce are also included in the same table. These a 



values are obtained by the Newton-Raphson procedure 
to zero the function /(a) = E(a) — E exact , where E(a) 
and E exact are the self-consistent atomic energy for each 
value of a and the exact total energies, respectively. The 
optimal a values are obtained for exact atomic energies 
in the highest symmetry for which the solutions have 
integral occupation numbers. The third set of a val- 
ues is obtained from a spherically-symmetric, numerical 
electronic-structure code for atoms. The a values for the 
DZVP2/A2 basis are usually larger than those for the 6- 
311G**/RI-J as the orbital variational principle requires 
lower energy from a larger basis set and increasing a low- 
ers the total energy. The fitting basis has a smaller effect 
on the total energies than the orbital basisA 

The second step consists of computing total energy 
of molecules using Eq. 1 and these a values, which in 
effect extrapolates the atomic energies. All molecules 
belonging to the G2 set are optimized and the abso- 
lute errors in the total energies of these molecules are 
computed. These are given in Table [H] We choose 
the G2 set as accurate experimental atomization ener- 
gies and thus total energies for these molecules are well 
known. We used the zero-point corrected atomization 
energies from Ref.0 and the atomic energies from Ta- 
ble I to synthesize an experimental total energy. The G2 
set is routinely used for performance appraisal of density 
functional models. While our energy is not obviously a 
density- functional, it shares a common roofei with DFT. 
Of all methods with this common root, only one other 5 
has been treated analytically in these fifty years. The 
errors in the total energies computed in our model are 
compared with the popular density functional and hybrid 
models in Table IHII These density functional models are 
certainly more sophisticated intrinsically than either the 
Gaspar-Kohn-Sham (GKS) or Xa methods. The total 
energies corresponding to these sophisticated functionals, 
however, are not obtained in particularly sophisticated 
fashion. It is necessary to use a numerical grid for these 
more sophisticated functionals. Consequently the ener- 
gies are dependent on how one (here GAUSSIAN03) 22 
chooses to orient each molecule relative to that neces- 
sary grid. On the other hand, energies obtained in our 
model use a less sophisticated functional, but are in- 
dependent of molecular orientation, and thus could be 
obtained to whatever accuracy is needed to judge be- 
tween analytic DFT models. The sophisticated pure and 
hybrid density functional models that we compare with 
are: the parameter-free PBEiS and the three-parameter 
hybrid density functional B3LYP, which mixes Hartree- 
Fock (HF) exchange with the local exchange functional 
and Becke's generalized gradient exchange functional 11 
along with local correlation^ and the LYP— correla- 
tion functionals. The mixing coefficients were empiri- 
cally determined to minimize atomization energies, not 
total (less zero-point) energies addressed in this work. 
We also computed the total energies for several other 
density functional models. These include the local den- 
sity approximation (LDA)ji£ BLYP-GGAfiiJ^ Perdew- 
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0.74324 
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0.74262 
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TABLE II: Absolute Error in total molecular energies with 
respect to its exact value for the G2 set of molecules. The 
computational models are- Ml: present/6311G**/RI-J, M2: 
PBE/6311G** M3: B3LYP/6311G*. The absolute errors are 
at the optimized geomertries in the repsective model. Last 
column contains "exact" total energy. The errors and mean 
absolute error (MAE) are in kcal/mol. (See text for more 
details) . 





Ml 


M2 


M3 


Exact (a.u.) 


H 2 


1.83 


0.44 


0.20 


-1.17530 


LiH 


1.44 


1.17 


0.75 


-8.06995 


BeH 


0.64 


1.43 


0.88 


-15.24604 


CH 


1.23 


2.48 


0.63 


-38.47838 


CH 2 ( 3 Si) 
CH 2 (Ui) 


0.40 


2.20 


0.78 


-39.14747 


1.64 


2.75 


0.54 


-39.13265 


CH 3 


0.53 


2.32 


0.94 


-39.83328 


CH 4 


0.87 


2.46 


0.96 


-40.51288 


NH 


1.20 


2.76 


0.65 


-55.22291 


NH 2 


1.79 


2.78 


0.73 


-55.87924 


NH 3 


1.87 


3.03 


0.57 


-56.56378 


OH 


0.56 


3.19 


0.84 


-75.73640 


H 2 


0.30 


3.53 


0.49 


-76.43687 


HF 


0.30 


3.96 


0.67 


-100.45522 


Li 2 


1.31 


1.96 


0.96 


-14.99488 


LiF 


0.65 


4.65 


1.30 


-107.43019 


C 2 H 2 


1.29 


4.34 


0.87 


-77.33589 


C 2 H4 


0.67 


4.38 


1.28 


-78.58624 


C 2 He 


0.08 


4.56 


1.55 


-79.82257 


CN 


0.84 


4.13 


0.77 


-92.72025 


HCN 


0.10 


4.79 


0.60 


-93.43906 


CO 


1.77 


4.93 


0.98 


-113.32506 


HCO 


2.14 


4.49 


1.43 


-113.85550 


H 2 CO 


1.55 


4.85 


1.35 


-114.50705 


H3COH 


0.75 


5.36 


1.39 


-115.72729 


N 2 


0.99 


5.02 


0.54 


-109.54414 


N 2 H 4 


1.78 


5.34 


1.10 


-111.87768 


NO 


0.76 


4.85 


1.20 


-129.90066 


o 2 


2.73 


4.73 


1.80 


-150.32587 


H 2 2 


1.08 


5.79 


1.38 


-151.56204 


F 2 


2.06 


6.38 


1.83 


-199.52335 


C0 2 


4.97 


6.81 


1.96 


-188.59875 


SiH 2 (Mi) 


1.91 


7.13 


1.32 


-290.61127 


SiH 2 ( 3 J B 1 ) 


0.70 


6.76 


1.37 


-290.57828 


SiH 3 


2.35 


7.20 


1.38 


-291.23127 


SiH 4 


2.92 


7.51 


1.47 


-291.88218 


PH 2 


1.96 


7.76 


0.95 


-342.51350 


PH 3 


2.84 


8.13 


0.79 


-343.15565 


H 2 S 


1.37 


9.56 


0.39 


-399.43051 


HC1 


0.36 


11.31 


1.65 


-460.86924 


Na 2 


0.87 


9.52 


2.47 


-324.54677 


Si 2 


0.19 


13.29 


1.82 


-578.85904 


P 2 


1.70 


15.59 


0.33 


-682.72677 


S 2 


0.63 


18.37 


1.35 


-796.44191 


Cl 2 


0.48 


22.60 


4.00 


-920.49227 


NaCl 


0.66 


16.28 


0.72 


-622.61585 


SiO 


0.57 


10.02 


1.10 


-364.74170 


CS 


0.59 


11.42 


0.54 


-436.25782 


SO 


1.28 


12.18 


0.32 


-473.40636 


CIO 


1.12 


14.00 


1.39 


-535.36947 


C1F 


1.41 


14.74 


1.25 


-560.02885 


Si 2 H 6 


3.98 


14.51 


2.52 


-582.58381 


CH 3 C1 


0.50 


13.37 


0.99 


-500.17224 


H 3 CSH 


0.36 


11.65 


0.20 


-438.73830 


HOC1 


0.75 


14.24 


1.35 


-536.02883 


so 2 


1.68 


16.03 


1.38 


-548.68595 


MAE 


17.3 


100.8 


15.3 





TABLE III: Mean absolute error (MAE) and mean error (ME) 
in total energies of G2 set of 56 molecules for different mod- 
els. The errors are in kcal/mol and are at optimized geome- 
tries in respective model. Ml: present/6-311G**/RI-J, M2: 
present/DGDZVP2/A2 (see text for more details). 



Model 


MAE 


ME 


present (Ml) 


17 


-0.4 


present (M2) 


16 


-0.5 


LDA 


532 


20.0 


PBE-GGA 


101 


3.9 


BLYP-GGA 


13 


0.4 


PW91-GGA 


19 


3.8 


PBE-GGA (NRLMOL) 


87 


3.1 


HCTH407-GGA 


14 


0.3 


VSXC-meta GGA 


60 


-1.7 


B3LYP-hybrid GGA 


15 


0.3 


PBElPBE-hybrid GGA 


91 


3.9 



Wang (PW91) GGk^ the meta-GGA functional con- 
taining kinetic densities due to Voorhis and Scuseria 
(VSXC)fi^ the empirical GGA due to Handy and Co- 
hen (HCTH407),2a and the hybrid PBE1PBE (also called 
PBE0) functional^ The HCTH407 functional contains 
15 parameters that are fitted to a database of 407 prop- 
erties including total energies. The PBE1PBE model is a 
hybrid functional of PBE and 25% of HF exchange. The 
total energies used for comparison are computed with the 
6-311G** basis set by the GAUSSIAN03 code at fully 
optimized geometries within each respective model. The 
default fine mesh is used in all calculation. 

It is well known that Gaussian are not ideal for com- 
puting atomic energies, and therefore larger Gaussian ba- 
sis sets must be used. Thus we also calculated total ener- 
gies for the PBE models at the geometries optimized by 
the NRLMOL code^ using the NRLMOL basis^i The 
NRLMOL Gaussian basis set is optimized for the PBE 
density functional model and is much larger than the 6- 
31 1G** basis. These results are also included in Table 
IIIII The mean G2-set errors are also tabulated as a test 
for this, and potentially other, systematic errors. As the 
mean errors are small, Gaussian basis set incompleteness 
in computing atomic energies is not an important con- 
sideration. 

In this study of extrapolating atomic energies only the 
MAE in total energy for the density functional models 
in Table IIIII is of concern. The zero-point energy is not 
included in total energies as the analytic second deriva- 
tive calculation required for vibrational frequencies is not 
yet available in our model. Correcting for zero-point mo- 
tion can introduce a small error particularly for poly- 
atomic systems. The MAE in total energy obtained by 
our method is 16.2 kcal/mol for the DZVP2/A2 and 17.3 
kcal/mol for the 6-311G**/RI-J basis. The comparison 
of our MAE with that obtained from the most widely 
used pure and hybrid density-functional models is favor- 
able. This simple way of determining total energies of 
molecules by extrapolation is remarkably accurate. It is 
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also likely to be significantly more accurate than the ab 
initio methods that can treat the largest G2 molecules 
today because the HF energy bounds the exact energy 
from above. 

Our model gives total energies comparable to the 
B3LYP, which is widely accepted as the most accurate 
hybrid density-functional model, because it gives the best 
atomization energies. The B3LYP/6-311G** MAE in at- 
omization energy (for the G2 set) is about 4.5 kcal/mol. 
Using the larger 6-311+G(3df,2p) basis set the B3LYP 
MAE in atomization energy drops to 2.27 kcal/mol. As 
the error in the total atomic energies in the present model 
is zero by construction, its error in atomization energy is 
the same as its error in the total energy. For B3LYP, and 
perhaps other models, this essentially single experimen- 
tal quantity, may be computed in ways that can differ 
by almost an order of magnitude in accuracy. Perhaps, 
density-functional models must benefit from cancellation 
of errors in total energies of atoms and molecules to give 
atomization energies comparable to or better than ab ini- 
tio methods. Unexpectedly, these models provide quite 
accurate total energies as well as atomization energies. 

The mixing parameters in the hybrid B3LYP method 
are empirically obtained by minimizing the errors in at- 
omization energies. Like B3LYP, our method can also be 
parameterized for atomization energies by optimizing the 
MAE in atomization energies of the G2 set of molecules. 
Such an optimized parametrization gives MAE in atom- 
ization energy that is intermediate between the PBE- 



GGA and B3LYP. The calculations could be further im- 
proved in several ways. A larger basis also lowers the 
B3LYP error A larger basis might lower our MAE, but 
the basis-set effect is largely canceled by adjusting a to 
get exact atomic energies, and high angular-momentum 
functions are less important for HF-less functionals. Per- 
haps the a 's should be adjusted towards the GKS value 
for valence or open-shell electrons. That seems algorith- 
mically possible, but departs from Slater's original model 
for the binding of atoms in molecules and solids. 

In conclusion, we have shown that fairly accurate es- 
timates of the total energy of molecules can be obtained 
by extrapolation using experimental atomic energies and 
an analytic quantum mechanical model. The computa- 
tional model is illustrated by computing the geometry- 
optimized total electronic energy of the G2 set of 56 
molecules, without grids and thus with machine-precision 
matrix elements. It is shown to perform as well as or 
better than some widely used pure and hybrid density- 
functional models in computing total molecular elec- 
tronic energies simply by extrapolating atomic energies 
quantum mechanically, according to Slater's prescription. 
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of Defense's High Performance Computing Moderniza- 
tion Program, through the CHSSI Project MBD-5, sup- 
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